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I. INTRODUCTION 



To understand the origin of CP violation (CPV) is an important subject not only for 
exploring the basic symmetry of space-time and elementary particles but also for under- 
standing the evolution of our universe. It is well known that in the Standard Model(SM) 
of particle physics, CP violation is characterized by a single weak phase in Cabibbo- 
Kobayashi-Maskawa matrix|l|], which can provide a well explanation for the direct CP 
violation e'/e established in kaon decays^, and also direct CP violation[4] observed in 
B-meson decays^. Though the theory of the strong and electroweak(EW) interactions in 
SM has met with extraordinary success, it is widely believed that the SM can not be the 
final theory of particle physics, in particular because the Higgs sector of SM is not well 
understood yet and the CP phase in CKM matrix is not enough to understand the baryon 
and anti-baryon asymmetry in the universe. It was suggested that CP symmetry may be 
broken down spontaneously^. Many possible extensions of SM in Higgs sector have been 
proposed^. Other possible extensions of the SM have been explored, such as the Super 
Symmetric model(SUSY), little Higgs model and extra dimensions, which all make better 
the situation of the SM. But no single model is good enough to solve all the problems exist- 
ing in the SM and it is then worthwhile to consider all the possibilities beyond SM. As one 
of the simplest extensions of the SM, the so-called two-higgs-doublet model (2HDM) which 
introduces an extra Higgs doublet without imposing the ad hoc discrete symmetries has 
been investigated widely from various considerations 0, 
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Motivated solely from the origin of CP violation, a general two Higgs doublet model with 
spontaneous CP violation (Model III 2HDM) has been shown to provide one of the sim- 
plest and attractive models in understanding the origin and mechanism of CP violation 
at weak scale (l2l. 13| . In such a model, there exists more physical neutral and charged 
Higgs bosons and rich induced CP-violating sources from a single CP phase of vacuum. 
Of particular, the model III 2HDM allows flavor- changing neutral currents but suppressed 
by approximate U(l) flavor symmetry, which is different from the so-called model I and 
model II 2HDM in which a ad-hoc discrete symmetry [Zi symmetry) has been imposed 
to avoid the FCNC. 

It is known that the FCNC's concerning the first two generations are highly suppressed 
from low-energy experiments, and those involving the third generation is not as severely 
suppressed as the first two generations. So the model III 2HDM can be parameterized 
in a way to satisfy the current experimental constraints. The constraints on Model III 
2HDM from neutral meson mixings (K° — K°, D° — D°, B° — B°) [l8| and radiative decays 
of bottom quark [IE [20, 21] have been studied in details. In this note, we shall investigate 
the possible new effects of model III 2HDM on two-body charmless nonleptonic B decays 
B — ► h\h2 with h\,h2 being the charmless light hadrons. This is because those decays 
have triggered considerable theoretical interest in understanding SM. Also those decay 
channels are also thought to be sensitive and important in exploring new physics beyond 
the SM as they involve the so-called tree (current-current) b — ► (it, c) and/or B — ► (d,s) 
penguin amplitudes with both QCD and electroweak penguin transition participating. 
In the 2HDM, there are five Higgs particles including the H° Higgs in SM, these extra 
Higgs will mediate all the penguin transitions. As the couplings involving Higgs bosons 
and fermions have complex CP phases in the model III 2HDM, CP violation effects occur 
even in the simplest case that all the tree-level FCNC couplings are negligible. With the 
improvement of experimental precision, more and more direct CPV have been observed 
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and will be much precisely tested in the future experiments. 

The paper is organized as follows. In Sec II, we first describe the theoretical frame 
including a brief introduction of the two-Higgs-doublet model with spontaneous CP vi- 
olation, i.e., Model III 2HDM, and the effective Hamiltonian as well as the generalized 
factorization formula, which is our basic tool to estimate the branching ratios and CPV 
asymmetry of B meson decays. In Sec III, we make a detailed calculation with numer- 
ical results evaluated from a factorization ansatz which allows us to express the matrix 
elements < hih 2 \H e jf\B > as a product of two factors < hi\Ji\B >< h 2 \J2\B >, and 
make quantitative predictions. Our conclusions and discussions are presented in the last 
section. 



II. THEORETICAL FRAMEWORK 

A. Outline of the Two-Higgs-doublet Model 

One of the important developments of SM is the so-called Higgs mechanism, i.e., a 
spontaneous symmetry breaking mechanism by which the gauge bosons and fermions 
can get their masses. In the SM, a single Higgs doublet of SU(2) is sufficient to break 
the SU{2) L x U(l)y symmetry to U(l) em and generate mass to the gauge bosons and 
fermions. Nevertheless, the Higgs sector of the SM has not been experimentally tested 
although enormous efforts have been made. For the origin of CP violation, SM gives no 
explanation as there is only one single neutral Higgs in SM and its interaction coupling 
constants are fixed by the known parameters and the fermion masses. Many attempts 
have been made by both theoretists and experimentalists to explore the mechanisms of 
CP violation since the discovery of CP violation in 1964. Spontaneous CP violation 
requires at least two Higgs doublets. A consistent and simple model which provides a 
spontaneous CP violation mechanism was constructed completely in a general Two-Higgs- 



doublet model[12l.ll3l|. Such a model III 2HDM not only explains the origin of CP violation 
in the SM, but also induces rich new resources of CP violation. The new sources of CP 
violation can lead to some new phenomenological effects which are promising to be tested 
by the future B factory and LHCb. In this note, we will focus on the phenomenological 
applications of the model III 2HDM in the two-body charmless hadronic B — > PV decays 



The two complex Higgs doublets in the Model III 2HDM are expressed as|12l. Il3l . 114 

MM- 



The corresponding Higgs potential can simply be written in the following general form: 
V{<j>) = A 1 (0 t 1 1 -^) 2 + A 2 (^ 2 -^ 2 2 ) 2 

+ A3(0l01 - ^1 2 )(0202 - \v\) + A 4 [(<M0 2 )(0 2 0l)] 

+^A 5 (<M02 + 2 0i - v i v 2 cos 5) 2 + A 6 O|0 2 - 2 0i - v ± v 2 sin 5) 2 

+ [\ 7 (44>i - \v 2 ) 2 + A 8 (4</>2 - \vl) 2 ][4>\<p2 + 4>\<Pl - i>iv 2 cos<J], (2) 
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where \{i = 1,2, ...8) are all real parameters. If all Aj are non-negative, the minimum 
occurs at : 



< 



< 



>= vie 10 , 
" >= v 2 , 



(3) 



With v±,V2 are the vacuum expectation values of 0i,02 respectively, and S the relative 
phase of the vacuum. It is clear that in the above potential CP nonconservation can only 
occur through the vacuum with 5^0. Obviously, such a CP violation appears as an 
explicit one in the potential when X 6 ^ o[l~3|. 

After a unitary transformation, it is natural and convenient to use the following basis: 



H 2 



1 

i 

71 



V2G+ 

v + 0? + iG° 

V2H+ 
0° + iA° 



(4) 



with: 



ve 




i8 



(5) 



where v = yvf+~v% and is related to the W mass by M w = gv/2. Here H° plays the 
role of the Higgs boson in the standard model. are the charged scalar pair with 
H ± = sin f3(j)fe~ %& — cosf3(J)f, where tan (3 = v 2 /vi. And as for the neutral Higgs, (j>x,(/>2 
are not the neutral mass eigenstates but linear combinations of CP-even neutral Higgs 
boson mass eigenstates, H and ho' 



H 
h 



cos a + 4>2 sin a, 
-0? sin a + 0° cos a i 



(6) 



where a is the mixing angle and when a = 0, ((p®, 0°) are identical with (Ho, ho) 
simplicity, the mixing with the pseudoscalar A is not considered here. 
Let us consider a Yukawa Lagrangian of the following form: 



For 



u 



where <f>i(i = 1,2) are the two Higgs doublets, .0^2 = ir 2 0i 2 , Qi,L{Uj,n) with i = (1,2,3) 
are the left-handed isodoublet quarks (right-handed up-type quarks), D^r are the 



righthanded isosinglet down- type quarks, while and £ 2 ij (hJ = 1)2,3 are family 
index ) are generally the nondiagonal matrices of the Yukawa coupling. After diagonal- 
izing the mass matrix of quark fields, the Yukawa Lagrangian that related to the decays 
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we considered in this paper can be written as: 

£ Y = --^—{H° cosa - h° sm^iJJMuU + DM D D) 

_ fl'1 na + fc' n ,a F({P 1 (1 +7j + {g ,l (1 _ is))u 

+D(f\(l + 7s) + « D, i(l - 75))0] 

ylO 1 1 1 1 

' \U{e-{\ + 75) - £^(1 - 7b))U - D(£ D -(l + 75) - e Dt -(l - 75))^] 



v^ L vs 2 V ' 7 s 2 V ' " vs 2 V ' ' s 2 

-#WgW d ~(i + 75) - ^ f ~(i - 7 5 )P 

-i?-D[^V^ M ^(l - 75) - V^ M ^(1 + 75 )]t/, (8) 

where [/ represents the mass eigenstates of u, c, t quarks and D represents the mass 
eigenstates of d,s,b quarks, Vckm is the Cabibbo-Kobayashi-Maskawa matrix and ^ U,D 
are the FCNC couplings in the mass eigenstate, and they may be parameterized in terms 
of the quark mass: 



v.d s // ^7777777" 
V2M W 



t u — i u ■ Vckm, 

i D = Vckm • £ D , (9) 

The first two generations' FCNC are naturally suppressed by the small quark masses,but 
the third generation has more space to get FCNC contributions. In this paper, we just 
choose the ^ U,D to be diagonal = t^' D (i = s, c, b, t) , and neglect the first generation 
quarks' contributions. So the really leading contribution arises from the diagram with a 
top quark in the loop and the relevant couplings will be ^ S ' D and they are explicitly 
given by: 

£ = &v u , g = tfv th , (10) 

From the above parameterization, the free parameters in this model are A^ (i, j = s, c, t, b). 
Their values can be constrained through experiments. 

In the model III 2HDM with spontaneous CP violation, the induced CP violation can 
be classified into the following four types via their interactions (l2l. fill ]: i) from the CKM 
matrix; ii) from the charged Higgs couplings to the fermions Charged] iii) from the neutral 
Higgs couplings to the fermions ^ neu trai] iv) from the CP nonconservation Higgs potential 
V{4>) via mixings among scalars and pseudoscalar bosons. 

The model allows flavor-changing-neutral-currents(FCNC) at tree level and via loop 
effects due to exchanges of Higgs bosons. One of the most stringent tests is from the 
radiative decay of B mesons and also from the inclusive decay rate of b — > 57 which has 
the least hadronic uncertainties. Other constraints could come from the B — B mixing, 
Po, Rb and the neutron electric dipole moment etc. In this note, we shall consider possible 
new effects in charmless hadronic two body decays of bottom mesons. 
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B. Effective Hamiltonian and Wilson coefficients 



The effective Hamiltonian for charmless B decays with AB 
G, 



1 is: 



* p=u,c i=3,. ..,16 

+Cj 1 Qj 1 + C^gQsg + C 7 Q 7 + C 8 Q 8 ) + h.c. 



Ill 



The operators Qi,...io, Q77, Qsg can be found in [24J, of which the Q\ and Q2 are the 
current-current operators and Q3 — Q6 are QCD penguin operators. Qj y and Q$ g are, 
respectively, the magnetic penguin operators for b — > 57 and 6 — > sg. Here the mass of 
the external strange quark is neglected compared to the external bottom-quark mass. 

The additional new operators related to the neutral Higgs mediated processes(6 — > sqq 
)are(25l|: 



Qu = 


(sb) {S+P )^2{qq)(S-P), 


Qn = 


{s i b j ) ( s+p)^2{q j q l )(s-p), 


Ql3 = 


(sb) ( s +P )^2(qq)(S+P), 


Qu — 


(sib j )( S +p)^2(q j qi)(s+p), 


Q15 — 


(5(7^(1 + 75)6)^(9^(1 + 75)9), 


<5l6 = 


{-§^{1 + 75)6,0 J2(g-y u (i + 75)%), 



(12) 

where {q\q2)s±p = <Ti(l ±75)92, <? = u,d,s,c,b. The operators in Eq (fTTT) are ob- 
tained from the Qi by exchanging L <-> i?. As the primed operators's contributions are 
suppressed by m^/m^, we shall neglect their effects in our present considerations. The 
Wilson Coefficients C^i = 1, ...10 have been calculated at LO[22|, 23] and NLO 24] in 
SM and also at LO in 2HDM 26, 27 . Here we list their initial coefficient functions in the 
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2HDM 27, 28 



Ct{M w 
C 2 (M W 
C 3 (M W 
C 4 (M W 
C 5 (M W 
C 6 (M W 



6n 



C 7 (M W 
C 8 (M W 
C 9 (M W 
C 10 (M W 

and the LO CV 7 , C$ g are sufficient: 



11 a s {M w ) 
2 4vr ' 
i _llaJMv] 
6 4tt 



35 a 
18 47 : 



a: 



2 ^ {E (x t ) + £ m (y)} + ^[25o(^) + C (x t )}, 



w 



8tt 

a s (M w 



24vr 



{^(x t )+^ w (y)}, 
{£ (* 4 ) + £ m (y)}, 



87T 



{^)+^(y)} 



[4C (z t ) + A)(z t )], 





a 
6tt 




— [4C (x t ) + D Q (x t ) + — n 

sin c/ 



-(10S (^)-4C (x t ))], 



(13) 



CV 7 (Mvy) 
08s (-Ww> 



2 6 
L>(:r t ) £(?/) 



A tt | 2 + 5(y)|A tt A 66 |e 



\ tt \ 2 + E(y)\\ tt \ bb \e 



w 



(14) 



2 6 

where x t = mf/M^r, and y = vr? t jMn±2. The Inami-Lim functions A,B,D,E are 

known in SM and 2HDM0: 

For the new operators Q(n,i2...i6)j the corresponding Wilson coefficients Cj, z = 11, ...16 
at leading order have been calculated in 25|, [29( : 



C u (M w ) 
Ci 3 (M w ) 

C 12 (M W ) = 14 (M W ) = C 15 (M W ) = Ci 6 (M w ) = 0, 



An m r A* T 



(15) 



Here the explicit expression of Cq 1 , Cq 1 can be found in [29| . 

For the B — ► processes, the Wilson coefficient functions must run from the 
scale to the scale of 0{m b ). For C\ — C W) the NLO corrections should be included. While 
for Cs g and CV 7 , LO results are sufficient. The details for the running Wilson coefficients 
can be found in Ref . [24] . As for the neutral Higgs boson induced operators, the one loop 



anomalous dimension matrices can be divided into two distangled groups [25 
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RL 





On 


012 


On 


-16 





012 


-6 


2 



(16) 



7 



and 





013 


014 


015 


016 


013 


-16 





1/3 


-1 


014 


-6 


2 


-1/2 


-7/6 


Ol5 


16 


-48 


16/3 





Ol6 


-24 


-56 


6 


-38/3 



As no NLO Wilson coefficients Ci,i = 11, 12, ...16 are available, we may just use the LO 
Wilson coefficients for a numerical estimation. 



C. Generalized Factorization formula 



For our present purpose, we may use the generalized factorization method [30l. 131 



32l . |33j to evaluate the hadronic matrix elements. We know that in full theory, the 
leading order QCD corrections to the weak transition is of the form a s ln(M^/ — p 2 ) for 
massless quarks, where p is the off-shell momentum of external quark lines and depends 
on the system under consideration. We can choose a renormalization scale /x and separate 
ln(M^/ - p 2 ) = ln(M^//i 2 ) + ln(/i 2 / - p 2 ). The first part ln(M^//i 2 )is included in the 
Wilson coefficients c(/j) and summed over to all orders in a s using the renormalization 
group equation, while the second part is due to the matrix element evaluations and is 
small. It is related to the tree matrix element via: 

(0(jl)) = g(tl)(0)tree (18) 

with: 

g(p) ~ 1 + a s (/x)( 7 In + c) (19) 

—p z 

where the \x dependence of the matrix elements is approximately extracted out to the 
function g(fi), that is: 

(We//) = c(f,)g(fi)(0} tree = c eff (0} tree (20) 

the effective Wilson coefficients c e ^ should be in principle renormalization scale indepen- 
dent. Thus it is necessary to incorporate QCD and EW corrections to the operators: 

(Oi(/i)> = [I+ ^-m s {n) + ^MliiiOU (21) 

with 

cf f = [I + ^-rh s T (v,) + £m/(/i)]^(/i), (22) 
The perturbative QCD and EW corrections to the matrices m s and m e from the vertex 



and penguin diagrams can be found in [33|, |34j, [35 



8 



Using the following parameterization for decay constant and form factors: 
< 0\A^\P(q) > = ifpq^,<0\V^\V(p,e) >= f v m v e^ 

(23) 

we arrive at 

x bp,v V \{q 2 q 3 ) v _ A \ >< P^y^B > = 2f v m v F^ p {m 2 v ){e ■ p B ), 
x bv,p _< P \(q 2 q 3 ) v _ A \ >< V\(q x b)v- A \B > = 2f P m v A^ v (m 2 P )(e ■ p B ), (24) 

Using the Fierz Transformation, 

(V — A)(V + A) — > -2(5 -P)(S + P), 

(V-A)(V- A) -> (V- A)(V- A) (25) 

One can easily obtain all the Qi,...io tree level matrix elements |30l. 321. For the new 



operators Qn,...i6, the additional factorization formulas are [36 

<V(k,e*)\qa^q'\0> = -i^K - e'Mft, 

< P(p)\qa^k"q'\B(p B ) > = % - {q 2 (p + p B )» ~ (m 2 B - m P )q^f p , (26) 

m B + mp 

with k = p B — P and q = p B — p. fy and f p are the tensor decay constant of vector 
meson and the tensor form factor relevant to B —>■ P decays, e* is the polarization vector 
of vector meson. The hadronic matrix element is given by 

— 2f ± f p m 2 
< V(k,e*)\q'a^q\0 X P(p)\qa, u b\B(p B ) > = Jy Jt V (e*-p B ), (27) 

m B + mp 

The tree level matrix elements of Q(n,i2,...i6) can be factorized as (b — > s for example): 

2 

< PV\Q U \B > = a u - ■ ■ T < P|(?g)y-A|0 >< V\(sb) V - A \B >, 

(m b + m s )(m q + m y/ ) 

< PV\Q 12 \B > = ~ ai2 < P\@q) V -A\0 >< V\{sb) V -A\B > 

= \cli2 < V\tfq) v _ A \Q >< P\{sb) v _ A \B >, 

9 

< PV\Q 1S \B > = -a 13 - ■ ^f— < P\(q'q) v -A\0 >< V\{sb) v _ A \B >, 

(m b + m s )(m q + m q <) 

j Tfl^ 

< PV\Q 14 \B > = --Oi4 7 ; rr~ ' r < P\(q'q) v -A\0 >< V\(sb) v _ A \B > 

2 (m b + m s )(m q + m q/ ) 

= ~ < \VY^ u q\0 >< P\sa^b\B >, 

< PV\Q 15 \B > = 2a 15 < \VYa" u q\0 >< P\sa^b\B >, 

< PV\Q ie \B > = -a 16 < \VY^s\0 >< P\qa^b\B > 

2 

= 6a w - ■ ^ < P\@s) v ^ A \0 >< V\{qb) v _ A \B >, 

(m b + m q )(m q , + m s ) 

(28) 
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with: 



an 


= en + ^, 


012 


= on + w , 

c 




ai3 


= C13 + W , 

c 


Ol4 


, C13 

= C14 + ^7, 

c 




015 


, c 16 
= C15 + jf, 

c 


0l6 


. c 15 

- cie + jf, 


(29) 



iVg is the effective color number relative to the new six operators,which is set to be 
universal in all the decay channels. In this paper we fix it to be N' c = 3 to estimate the 
neutron higgs effects. As for the SM operators, besides the perturbative QCD and EW 
corrections to the hadronic matrix elements that can be factorized into the effective Wilson 
coefficients, there still exists the nonfactorizable effects, such as the spectator quark effects, 
annihilation diagrams and space-like penguins. Consider an arbitrary operator of the form 
O = qi a Tq2 q^^Vq" which arises from the Fierz transformation of a singlet- singlet operator 
with T and T' being some combinations of Dirac matrices. By using the identity: 

O = yiTq 2 q 3 T'q 4 + -fcA'T©^^^ (30) 
the matrix element of M — > P\P 2 can be expanded as: 

< P x P 2 \0\M > = \ < n|girg 2 |0 X P 2 \q 3 T'q A \M> f +^ < P^Tq^O >< P 2 \q 3 T'q 4 \M > nf 
+ ^ < P 1 P 2 \q 1 X a Tq 2 q 3 X a r'q4M > . (31) 

The last two terms on r.h.s are nonfactorizable, and their contributions are included in the 
effective color number N^f . To evaluate the decay amplitudes, it is useful to introduce 
the combination of Wilson coefficients 

21 ~ 2i + [m")**- 1 ' 



"2i-l - ^2i-l -r ^ f j 

The values of N^^ can be found in [321, that is: 



e ff _ e// I r e ff (QO\ 

)2i-X 



N ejf(y — A) = (iV^)i « (K ff h « (K") 3 « (N e J% « (Nff) 9 « (N; 



110, 



N ff(y + A) = (N e /f) 5 « (N e /f) 6 « {N*Jf) 7 « (N e /f) 8 , (33) 

As shown in [32| that in general N^^(V — A) 7^ N^^(V + A). The satisfied choice 
is that Nf f (V - A) < 3 < Nf s (V + A). And it is reasonable to take the value of 
N%ff(V — A) = 2,N^f(V + A) = 5. From now on, we will drop the superscript "eff" 
through the paper for convenience. 



III. B -» PV DECAYS IN MODEL III 2HDM 

Based on the effective Hamiltonian obtained via the operator product expansion and 
renormalization group evaluation, one can write down the amplitude for B — > PV decays 
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TABLE I: Input parameters 









M B , 


nib 


1.528 x 10 -12 pfl 


1.472 x 10~ l2 ps 


b.2&GeV 


5.37GeV 


4.2GeV 


mt 


m u 


m d 


m c 


m s 


!7AGeV 


3.2MeV 


QAMeV 


l.lGeV 


0A05GeV 


7T 




m„ 


m„' 


p 


O.UGeV 


0A35GeV 


0.5A7GeV 


0.958Gey 


0.77GeV 


m ± 




m0 


m K ± 


m K o 


0.77GeV 


0.782GeF 


1.02GeF 


OAMGeV 


0A98GeV 


m K *± 




Aqcd 


U 


Ik 


0.892Gey 


0.896Gey 


225MeF 


0.132Gey 


0A6GeV 


f P 




fx* 


u 


f T 

Jp 


0.2lGeV 


0.195Gey 


0.221GeF 


0.237GeV 


0AA7GeV 




Ik* 


J 4, 






0.133GeV 


0A56GeV 


0.183GeF 







and calculate the branching ratios and CP violating asymmetries once a method is derived 
for computing the hadronic matrix elements. For purpose of this paper, we are going to 
explore the new physics contributions to the exclusive decays B — > PV in the general 
model III 2HDM with spontaneous CP violation. For a numerical estimation, we will 
employ the generalized factorization approach described in the previous section. 

We begin with the following definitions for the branching ratio and CP violation asym- 
metry: 



-4 



CP 



\A\ 2 -\A\ 2 
\A\ 2 + \A\ 2 ' 

„3 



BR{B -> PV) = \^~rr B {\A\ 2 + \A\ 2 )/(e ■ p B f, (3 I ) 



v 



where A and A are the decay amplitudes of B and B respectively, e is the polarization 
vector of the vector meson. The input parameters in calculation are listed in Table.I. 

Here ju and fj^ are all decay constants of the mesons, Jm comes from the experimental 
measurements, but fL is calculated from quenched lattice QCD and QCD sum rulesjl^ 



38|.As for the form factor of pseudoscalar and vector mesons, we use the results from 



Light Cone Sum Rules (LCSR)[36|, [39|, but for the form factor of 7/, we use the value 



of BSW model[40(. And for the 77 — 77' mixing effects, we use the results of [41]. The 
B — > P(V) form factors' values are listed in Table. II. For comparison, we list both the 
results for light-cone sum rules(LCSR) and from sum rules in the framework of heavy 
quark effective field theory j42|. 

In the model III 2HDM, Xij(i,j = c, s,b,t),m H ±,m ho ,mA ,m Ho are free parameters 
that should be constrained from experiments. It was shown from B® s — B® s mixing that 
the parameters |A CC | and |A SS | can reach to be around lOOjllJ, and their phases are not 
constrained too much. In our present considerations, we simply fix their phases to be tt/4 
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TABLE II: The relevant Form Factors at q 2 = for B —* P transitions from LCSR 36j, 1371 . 
39] (the first row), sum rule in Heavy Quark Effective Field Theory 42] (the second row)and BSW 



model 



40j(the third row). The values in the square brackets are the B — > rj' form factors. 





Decay Channel 


B q -tt 


B q ^K 


B q - rP 


B S ^K 


B s - /> 




LCSR 


0.258 


0.331 


0.275[-] 








SRHQEFT 


0.285 


0.345 


0.247[-] 


0.296 


0.281[-] 




BSW 


0.333 


0.379 


0.307[0.254] 


0.274 


0.335[0.282] 



TABLE III: The relevant Form Factors at q 2 = for B — ► V transitions from LCSR[36J, 137 
[i^] (the first row), sum rule in Heavy Quark Effective Field Theory |42| (the second row)and BSW 
model[40](the third row). 





Decay Channel 


Bg^P 


B q -> u 


Bq — > K* 


B s ^4> 


B S ^K* 




LCSR 


0.303 


0.281 


0.374 


0.474 


0.363 




SRHQEFT 


0.363 


0.341 


0.400 


0.397 


0.337 




BSW 


0.281 


0.280 


0.321 


0.475 


0.364 



to see their effects. For X u and A&&, the constraints come from the experimental results of 
B — B mixing, T(b — > S7),rfb — » ctv t \ Pn, Rh and the electric dipole moments (EDMS) 
of the electron and neutron [3, 0, HH, 29, 43]. For a numerical calculation, we are going 
to consider the following three typical parameter spaces which are allowed by the present 
experiments: 



Case 


A : 


|A tt | 


= 0.15; 


A&6 


= 50, 


Case 


B : 


|A«| 


= 0.3; 


Abb 


= 30, 


Case 


C : 


|A«| 


= 0.03; 


\Kb 


= 100, 



and: 



a + bb = vr/2, 



For the Higgs mass, the following values are assumed: 

m Ao ^ 120GeV, m ho ~ 115GeV, 
m Ho ~ 160GeV, m H ± ~ 200GeV 

All the numerical results are presented in Table V ~ Table IX. 



(35) 



(36) 



(37) 



IV. CONCLUSIONS AND DISCUSSIONS 

As the charged Higgs mediated one loop FCNC effects to the AB = 1 charmless 
decays are mostly characterized through the Wilson coefficient , which is included 
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in the C^ i 5 6 7 and there are no new operators beyond the basic operators Qi,2...i0) 
their contributions to the Wilson coefficients are given in Table. III. On the contrary, 
the neutral Higgs mediated processes will bring in new operators Q(ii,i2,-,i6) with the 
new Wilson coefficients C(n,i2,---,i6)- They are nonzero when the neutral Higgs couples 
to the second and third generation of quarks, and the numerical results are presented in 
Table. IV. From the above calculations, it is seen that in some decay channels, the new 
physics contributions can be significant, especially to CP violations. 

a). As we have set the Yukawa couplings X iu and to be zero, so the neutral Higgs 
contributions to B — > (p,uj, K*)ir, (p,u)K decays are actually ignored, only the charged 
Higgs give new contributions. One can see that the branching ratio of B — > K°p° decay in 
the model III 2HDM is the same as SM prediction (about 1.55 x 10~ 6 ), which is far below 
the large central value of experimental result (5.4 ±0.9) x 10~ 6 . Though the annihilation 
diagram and exchange diagram are not taken into account, their contributions are still 
not enough to give such an enhancement. So one needs to find some new mechanism 
to explain this discrepancy. The same situation also appears in the B — > K~ p + , K*~tt + 
decays, where the experimental results are much larger than theoretical predictions both 
in SM and 2HDM when simply using the generalized factorization approach. Though the 
branching ratios could be enhanced by using improved QCD Factorization(QCDF) jlij . 
the resulting values are still smaller than the measured results, b). The model III 2HDM 
prediction for the CP violation of Bj — > K<p decay is 5 ~ 7 times larger than the SM 
prediction, which can be a signal to look for new physics in future experiments. But the 
prediction of the branching ratio are both smaller than the experimental one. c). The 
SM and model III 2HDM predictions for branching ratio of Bd — > K*°n° are the same in 
size and all consistent with the experimental result at la level. While the new physics 
prediction for CPV can flip the sign of the SM one and be 1 ~ 5 times larger in size 
and still within la error of the experiments, d). In B^ — > K*(j],T]') decays, new physics 
effects to CPV becomes significant. In B — > K*r], the 2HDM prediction is negative but 
the SM one is positive. In B — > K*r]', the 2HDM prediction can be as large as 40%, which 
is about seven times of SM prediction, e). In Bd — > p + n~ decay, the model III 2HDM 
prediction can enhance the CP violation from about —20% in SM to about —30%. Both 
SM and 2HDM predictions for the branching ratio of Bd — > p°7r° decay are much smaller 
than the experimental result. Such an inconsistence cannot be improved even in QCD 
factorization method[i3]. As for CP violation, the SM and model III 2HDM predictions 
have opposite sign with the magnitude (10 ~ 15)%. As the current experimental error is 
still too big to draw a conclusion, much more precise measurement is needed to test it. 
f). In Bd — > urn decays, the new physics effects to CP violation may be distinct with the 
SM prediction, as it not only flips the sign but also enhances the magnitude by a factor 
of three. 

For B® — > PV decays, the new physics contribution can be large in some decay chan- 
nels, a). In B s — > K*i], the 2HDM prediction enhanced the direct CP violation to about 
—50% compared to the SM prediction —28.8%, but for B s — > K*r)', the new physics 
contribution is destructive and reduce the SM prediction —37% to about —20%. b). In 
B s — > prp ) decays, new physics contribution to branching ratio is destructive but gives 
an enhancement to the CPV to about four times of the SM prediction, c). In B s — > ) 
decays, new physics effects to branching ratios and CPV are both significant, d). In 
B s — > K°K°* decay, 2HDM can give about 25 ~ 70% enhancement to the branching 
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ratio, e). In B — > K°(f) decay, new physics effect to CPV is very significant, the SM 
prediction is almost zero, but the new physics effects can enhance it to about —10%. 

For B u — > PV decays, there are also some new effects from the extra Higgs contribu- 
tions: a). In B u -> n-R* decay, the 2HDM prediction of CPV can be ten times of the 
SM one and are more closer to the experimental value, b). In B u — > K~cf) decay, CPV 
can be about 10% in 2HDM , which is much larger than the SM prediction 1.44% and is 
within 2a level of experimental results, c). In B u — > K*~7] decays, new physics contribu- 
tion reduced the CPV to about a half or a quarter of the SM one and is much closer to the 
experimental result, d). In B u — > p~r] decay, 2HDM prediction for CPV is 2 ~ 3 times of 
the SM one and much closer to the experimental central value, e). In B u — > ir~(j) decay, 
new physics enhancement for the branching ratio and direct CPV can be all significant, 
f). In B u — > K*~K° decay, 2HDM predictions for CPV are 20 ~ 24%, which is very much 
larger than the SM prediction —1.73% On the contrary, in B u — > K*°K~ decay, 2HDM 
prediction for CPV can be much smaller than the SM one. 

From the above results, we see that in some decay channels, the theoretical predictions 
for branching ratios are still far from the experimental results not only in SM but also 
in model III 2HDM,such as B —>■ Kp, K*tc decays. And even using the improved QCDF, 
the situation cannot be improved much. There must be some new mechanism to improve 
those situations. For simplicity, we have not considered the possible effects of final state 
interaction (FSI) and the contributions from annihilation and exchange diagrams although 
they may play a significant rule in some decay channels. As for factorization part, in 
principle, N%** can vary from channel to channel as in the case of charm decay. However, 
in the energetic two-body B decays, is expected to be process insensitive |30l. |32| . and 
the preferred values are obtained from the data to be N*ff(V — A) = 2, N%ff(V + A) = 
In a numerical calculation, we have considered only three cases for parameter 
choice in a general model III 2HDM to be consistent with the experimental results. Also we 
have totally neglected the first generation Yukawa couplings and the off-diagonal matrix 
elements of the Yukawa coupling matrix, such as Xtc,sb etc. to eliminate the FCNC at tree 
level. However, it is still possible that FCNC involving the third generation quarks exists 
at tree level, so the constraints can be less stronger to get nonzero off-diagonal elements. 

In conclusion, we have shown that the new Higgs bosons in the general model III 2HDM 
with spontaneous CP violation can bring out some significant effects in some charmless 
B-meson decays, which can be good signals in the future B factory experiments to test 
the SM and look for new physics from more precise measurements. 
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APPENDIX A: THE EFFECTIVE WILSON COEFFICIENTS 



TABLE IV: The effective Wilson coefficients C e ^[ 10 ) in b ~> s process in SM and 2HDM at 
i x = m h = A.2GeV 



Model 


SM 


Case A 


Case B 


Case C 


cf 


1.17 


1.17 


1.17 


1.17 




-0.37 


-0.37 


-0.37 


-0.37 




0.024 + 0.0035/ 


0.024 + 0.0067 


0.024 + 0.00487 


0.024 + 0.00757 


cf 


-0.050 - 0.0107 


-0.05 - 0.0187 


-0.05 - 0.0147 


-0.05 - 0.0237 


cf 


0.015 + 0.00357 


0.015 + 0.0067 


0.015 + 0.0057 


0.015 + 0.00757 


r eff 


-0.064 - 0.0107 


-0.064- 0.0187 


-0.064- 0.0147 


-0.064 - 0.0237 


cf 


-0.00028 - 0.000247 


-0.00035 - 0.000247 


-0.00035 - 0.000247 


-0.00035 - 0.000247 


r eff 


0.00055 


0.00061 


0.00061 


0.0006 


r eff 


-0.011 - 0.000247 


-0.011 - 0.000247 


-0.011 - 0.000247 


-0.011 - 0.000247 


r eff 
°10 


0.0038 


0.0034 


0.0034 


0.0034 



TABLE V: The Wilson coefficients Cf{{ 12 lfn at fi = m b = A.2GeV 



Parameter Space 


Case A 


Case B 


Case C 




-0.089 + 0.127 


-0.089 + 0.197 


-0.11 + 0.137 













Cl3 


-0.031 - 0.0517 


-0.054 - 0.0727 


-0.030 - 0.0557 




-0.00063 - 0.00107 


-0.0011 - 0.00157 


-0.000061 - 0.00117 


°15 


0.00035 + 0.000577 


0.00061 + 0.000807 


0.00034 + 0.000627 


Cl6 


-0.0011 - 0.001757 


-0.0019 - 0.00257 


-0.0010 - 0.00197 


(~<S 

°11 


-0.0085 + 0.0127 


-0.0085 + 0.0187 


-0.010 + 0.0127 













C*13 


-0.0030 - 0.00497 


-0.0052 - 0.00697 


-0.0029 - 0.00527 




-0.000060 - 0.000107 


-0.00011 - 0.000147 


-0.000059 - 0.000107 




0.000033 + 0.0000557 


0.000058 + 0.0000787 


0.000032 + 0.0000597 


Cl6 


-0.00010 - 0.000177 


-0.00018 - 0.000247 


-0.0001 - 0.000187 
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APPENDIX B: NUMERICAL RESULTS OF B — > PV DECAYS 



TABLE VI: CP averaged branching ratios (in units of 10 6 ) (first line) and direct CPV (second 
line)for charmless B° d -► PV decays in SM and 2HDM. Nt ff (V-A), Nt ff (V-A) are fixed to be 
2 and 5 respectively and N' c = 3. The parameter spaces are: Case A: (\Xtt\ = 0.15, |Abb| = 50, 6 = 
it/2), Case B: (|A tt | = 0.03, |A»| = 100,0 = tt/2), Case C: (|A«| = 0.3, |A»| = 30,6 = ir/2). And 



Decay channel 


SM 


Case A (2HDM) 


Case B (2HDM) 


Case C (2HDM) 


Exp 




1.56 


1.55 


1.55 


1.56 


5.4 ±0.9 




1.7% 


2.10% 


1.97% 


2.20% 


- 




2.01 


1.94 


1.97 


1.91 






-3.6% 


-3.83% 


-3.90% 


-3.76% 


(17^)% 


5 o _^ #-*-„.+ 


3.24 


3.92 


3.56 


4.33 


9.8 ± 1.1 




24.5% 


27.5% 


26.2% 


28.2% 


(-5 ± 14)% 


#o ^ K *o n o 


1.47 


1.47 


1.46 


1.51 


1.7 ±0.8 




-2.50% 


7.20% 


2.30% 


11.5% 




B° d - i^V 


4.80 


5.22 


5.23 


5.18 


».3±i:g 




1.40% 


5.96% 


10.0% 


10.3% 




#0 _^ K * v 


9.41 


10.4 


10.7 


10.8 


16.1 ± 1.0 




1.86% 


-2.68% 


-3.85% 


-1.93% 


(19 ± 5)% 


BjJ -► i^*??' 


1.33 


1.18 


1.49 


1.20 


3.8 ± 1.2 




5.50% 


32.7% 


22.1% 


40.4% 


(-8 ± 25)% 


£° - if°u; 


0.43 


0.44 


0.44 


0.44 


4.8 ±0.6 




0.00 


0.33% 


0.32% 


0.32% 




£° _> 


15.8 


15.3 


15.1 


15.1 


24.0 ±2.5 




-4.3% 


-4.4% 


-4.3% 


-4.4% 
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TABLE VII: Continue Table III 



Decay channel 


SM 


Case A (2HDM) 


Case B (2HDM) 


Case C (2HDM) 


Exp 


B° d - p+TT- 


17.3 


16.3 


17.7 


15.0 


24.0 ±2.5 




-18.8% 


-26.5% 


-26.4% 


-26.4% 




(I 


0.23 


0.24 


0.24 


0.25 


< 1.9 




-11.5% 


-10.1% 


-14.6% 


-6.45% 





B? -► K°K*° 

(I 


0.038 


0.037 


0.036 


0.037 


_ 




-1.73% 


20.3% 


22.0% 


24.3% 


_ 




0.0039 


0.0038 


0.0038 


0.0038 


< 0.6 




1.13% 


1.35% 


1.25% 


1.45% 




B° d </>r/ 


0.0023 


0.0022 


0.0024 


0.0022 


< 1.0 




1.13% 


1.35% 


1.25% 


1.45% 







0.015 


0.0017 


0.0017 


0.0016 


< 0.28 




3.90% 


1.35% 


1.25% 


1.45% 






0.80 


0.81 


0.82 


0.80 






-10.5% 


14.4% 


13.7% 


14.9% 


(-49+™)% 


B^PV 

(1 1 1 


0.82 


0.88 


0.83 


0.92 


< 1.5 




12.3% 


6.93% 


3.21% 


6.91% 





B° d w' 


0.50 


0.55 


0.54 


0.57 


< 3.7 




5.88% 


6.33% 


6.87% 


6.20% 




B° - u;tt 


0.60 


0.56 


0.53 


0.59 


< 1.2 




-4.97% 


12.4% 


12.7% 


12.1% 




S° -► wtj 


0.84 


0.77 


0.80 


0.73 


< 1.9 




-13.9% 


-11.1% 


-9.04% 


-11.2% 




5° -► wt/ 


0.53 


0.46 


0.50 


0.44 


< 2.8 




-19.7% 


-25.0% 


-26.0% 


-26.1% 
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TABLE VIII: CP averaged branching ratios (in units of 10 6 ) (first line) and direct CPV (second 
line)for charmless B® — > PV decays in SM and 2HDM. 



Decay channel 


SM 


Case A (2HDM) 


Case B (2HDM) 


Case C (2HDM) 


B Q S -► K* + 7r~ 


8.34 


8.73 


8.34 


8.44 




-0.13% 


-0.12% 


-0.13% 


-0.13% 


B° s - K+p- 


26.7 


27.7 


26.5 


26.2 




-4.3% 


-4.3% 


-4.3% 


-4.4% 


B ° _> K*°7T 


0.21 


0.20 


0.21 


0.21 




5.0% 


5.1% 


5.1% 


5.1% 


B° s - 


0.59 


0.70 


0.66 


0.73 




16.5% 


14.2% 


14.7% 


13.6% 


B° - 


0.70 


0.73 


0.69 


0.68 




-18.3% 


-18.3% 


-18.4% 


-18.1% 


5° -► K*?? 


0.29 


0.32 


0.31 


0.33 




-28.8% 


-43.2% 


-49.9% 


-42.8% 


5° -► X"V 


0.20 


0.23 


0.21 


0.25 




-37.1% 


-21.2% 


-21.3% 


-17.7% 


£° -► K~K* + 


1.98 


2.27 


2.24 


2.31 




-3.6% 


-3.3% 


-3.4% 


-3.2% 


5° -► i-C+iT^ 


5.49 


6.97 


6.81 


7.02 




24.5% 


22.8% 


20.3% 


21.6% 


B° - OT 


0.21 


0.04 


0.04 


0.04 




3.6% 


18.4% 


18.4% 


18.3% 


B° - 


0.12 


0.03 


0.02 


0.03 




3.6% 


18.3% 


18.4% 


18.4% 


5° - 07T 


0.21 


0.04 


0.03 


0.04 




0.00% 


0.00 


0.00 


0.00 


B° - wtj 


0.02 


0.02 


0.02 


0.02 




43.1% 


40.0% 


39.6% 


40.1% 


5° -► LOT]' 


0.01 


0.01 


0.01 


0.01 




43.1% 


40.0% 


39.6% 


40.1% 


B° s - <^ 


12.2 


17.9 


17.7 


20.0 




2.2% 


-12.3% 


-10.6% 


-14.2% 




0.61 


1.68 


1.92 


2.23 




15.6% 


-21.2% 


-21.3% 


-17.7% 


5° -► 


5.78 


8.03 


7.00 


9.14 




1.27% 


4.10% 


2.29% 


5.0% 


5° -► 


1.28 


1.03 


1.00 


1.06 




1.0% 


1.2% 


2.2% 


1.0% 


B° <pK° 


0.14 


0.12 


0.12 


0.12 




0.23% 


-9.7% 


-8.8% 


-11.7% 
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TABLE IX: CP averaged branching ratios (in units of 10 6 ) (first line) and direct CPV (second 
line)for charmless B~ — > PV decays in SM and 2HDM. 



Decay channel 


SM 


Case A 2HDM) 


Case B (2HDM) 


Case C (2HDM) 


Exp. 


B~ - K-p° 


0.59 


0.58 


0.58 


0.56 


4.25±8:™ 




-1.88% 


-2.40% 


-2.32% 


-2.38% 


(31115)% 


B- - K*~tt° 


2.62 


2.80 


3.25 


2.80 


6.9 ±2.3 




18.9% 


22.3% 


20.4% 


22.3% 


(4 ±29)% 


B~ - KOp- 


1.09 


1.10 


1.10 


1.10 


< 48 




0.34% 


1.15% 


0.75% 


1.15% 


— 




3.67 


3.88 


3.66 


3.88 


11.3 ±1.0 




-1.28% 


-12.9% 


-5.41% 


-12.9% 


(-8.6 ±5.6)% 


B" - K-o; 


2.22 


2.17 


2.20 


2.17 


6.9 ±0.5 




0.00% 


-0.33% 


-0.32% 


-0.33% 


(5 ±6)% 


B" - K~<f> 


5.16 


5.92 


5.57 


5.93 


8.30 ± 0.65 




1.44% 


11.7% 


10.0% 


10.3% 


(3.4 ± 4.4)% 


B~ -+ K*- n 


9.36 


10.51 


11.09 


10.6 


19.5±\£ 




13.6% 


6.73% 


10.2% 


4.52% 


(2 ± 6)% 


B- -+ K*~rf 


1.53 


1.32 


1.38 


1.37 






52.2% 


55.4% 


52.4% 


56.9% 


(30l? 7 )% 


B~ - 7T°p- 


11.4 


11.1 


11.3 


11.1 


ICS!};* 




-3.0% 


-3.1% 


-3.0% 


-3.1% 


(2 ±11)% 


B~ - vr-p 


7.36 


7.75 


7.50 


7.75 






4.2% 


4.1% 


4.2% 


4.1% 


(-7t\l)% 


B~ -► 


6.85 


6.50 


6.65 


6.50 


6.7 ±0.6 




-4.7% 


-4.8% 


-4.7% 


-4.8% 


(-4 ± 7)% 




11.1 


13.2 


11.0 


10.8 


5.3+_\i 




-0.90% 


-2.35% 


-3.34% 


-2.36% 


(1 ± 16)% 


- 


14.0 


12.7 


13.7 


12.8 


q 1 +3.7 




-9.9% 


-10.1% 


-9.60% 


-10.1% 


(-4 ± 28)% 


B" - 


0.0036 


0.016 


0.016 


0.016 


< 0.24 




1.13% 


15.5% 


7.96% 


15.5% 




B~ -»• iT^K 


0.038 


0.037 


0.036 


0.037 






-1.73% 


20.3% 


22.0% 


24.3% 




B~ -► if-K* 


0.25 


0.27 


0.26 


0.27 


< 5.3 




-37.1% 


-5.13% 


-14.6% 


-6.45% 
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